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ABdThACT 


'Satellite  measurements  have  now  confirmed  that  intense  radio  emis¬ 
sions  called  auroral  kilometric  radiation  are  generated  at  altitudes  of 
2  to  3  over  the  auroral  regions.  These  radio  emissions  are  believed 
to  be  generated  by  Doppler-shifted  cyclotron  radiation  from  the  elec¬ 
trons  involved  in  the  auroral  acceleration  process.  Using  current 
theories  for  the  generation  of  this  radiation  we  discuss  the  possibil¬ 
ities  for  obtaining  information  on  the  auroral  acceleration  processes 
from  the  spectrum  of  this  radiation.  For  example,  under  certain  condi¬ 
tions  it  is  shown  that  the  low  frequency  cutoff  provides  a  direct  indi¬ 
cation  of  the  upper  altitude  limit  of  the  acceleration  region.  Also, 
certain  drifting  features  in  the  spectrum  can  be  interpreted  in  terms  of 
the  propagation  of  shock-like  disturbances  along  the  auroral  field  lines 
at  velocities  near  the  ion-acoustic  speed. 


I.  INTRODUCTION 


Spacecraft  observations  over  the  last  decade  have  revealed  that  the 
earth  is  an  intense  radio  source  in  the  kilometer  wavelength  range  and 
that  the  generation  of  this  radiation  is  closely  associated  with  the 
auroral  acceleration  process.  The  first  reports  of  this  radiation  were 
by  Benediktov  et  al.  [1965I  and  Dunkel  et  al.  [ 197 0 ]  who  identified 
intense  radio  emissions  which  were  closely  correlated  with  high  latitude 
magnetic  disturbances.  Later  studies  by  Gurnett  [1971*]  showed  that  the 
peak  intensity  occurred  at  kilometer  wavelengths  and  that  the  radio  emis¬ 
sions  were  generated  at  low  altitudes  over  the  auroral  regions  in  asso¬ 
ciation  with  discrete  auroral  arcs.  During  bright  auroral  displays  the 
maximum  power  radiated  was  found  to  be  very  large,  ~  10  9  vatts.  The 
earth  is,  therefore,  a  very  intense  radio  emitter,  comparable  in  many 
respects  to  Jupiter,  which  has  long  been  recognized  as  an  intense  plane¬ 
tary  radio  source.  Because  of  the  association  with  discrete  auroral 
arcs,  Gurnett  [197*+1  suggested  that  the  radiation  was  associated  with  the 
inverted-V  electron  precipitation  bands  detected  by  low-altitude  satel¬ 
lites  [Frank  and  Ackerson,  1971]  •  This  relationship  has  now  been  con¬ 
firmed  by  direct  in  situ  measurements  in  the  auroral  regions  [Bensen  and 
Calvert,  1979i  Green  et  al.,  19791.  Polarization  measurements  by  Gurnett 
and  Green  [1978]  and  Kaiser  et  al.  [1978]  showed  that  kilometric  radi¬ 
ation  is  generated  in  the  right-hand  extraordinary  mode.  Of  the  many 
theories  which  have  been  proposed,  the  observed  polarization  and  the 
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absence  of  intense  electrostatic  waves  in  the  source  region  have  greatly 
restricted  the  possible  generation  mechanisms.  The  most  promising  mech¬ 
anisms  at  the  present  time  appear  to  be  the  direct  electromagnetic 
instabilities  proposed  by  Melrose  [197&1  and  Wu  and  Lee  [ 197 9 1  - 

In  this  paper  we  consider  the  relationship  of  the  auroral  kilo- 
metric  radiation  spectrum  to  the  spatial  structure  and  processes  occur¬ 
ring  in  the  auroral  acceleration  region.  The  basic  objective  is  to 
explore  the  possibilities  of  using  the  emission  spectrum  to  provide 
information  about  the  details  of  the  acceleration  process.  The  methods 
employed  may  also  have  much  broader  applications.  It  now  appears  that 
the  radio  emission  mechanism  operative  in  the  terrestrial  auroral 
regions  is  a  universal  process  occurring  in  association  with  the  forma¬ 
tion  of  auroral  arcs  throughout  the  solar  system.  For  example,  at  both 
Jupiter  and  Saturn  the  same  basic  radio  emission  process  is  believed  to 
be  associated  with  auroral  arcs  occurring  at  these  planets  [Warwick  et 
al.,  1979i  Kaiser,  19801  .  Furthermore,  radio  emissions  of  the  same  type 
have  also  been  tentatively  identified  from  the  planet  Uranus  [Brown, 
1976]  .  Since  it  will  not  be  possible  to  conduct  direct  in  situ  measure¬ 
ments  of  aurora  at  any  of  these  planets  in  the  near  future,  any  advances 
which  can  be  made  in  understanding  the  relationship  of  the  radio  emis¬ 
sion  spectrum  to  the  auroral  acceleration  may  provide  a  valuable  tool 
for  studying  auroral  acceleration  processes  in  other  regions  of  the 
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II.  UPPER  AND  LOWER  CUTOFF  FREQUENCIES 


A  basic  characteristic  (if  the  kilometric  radio  emission  is  the 
existence  of  well-defined  upper  .and  lower  cutofF  frequencies.  Figure  1, 
for  example,  shows  a  typical  auroral  kilometric  radiation  event  observed 
by  the  ISEE  1  spacecraft  in  the  local  evening  region  of  the  magneto¬ 
sphere.  As  can  be  seen  the  radio  emission  spectrum  has  clearly  identi¬ 
fied  upper  and  lower  cutoffs  which  fluctuate  over  a  wide  range  on  a  time 
scale  of  ten  minutes  or  less.  Except  for  fine  structure,  which  will  be 
discussed  later,  the  radio  emission  normally  extends  continuously  over 
the  entire  frequency  band  between  the  upper  and  lower  cutoffs.  Occasion¬ 
ally  the  radiation  disappears  for  short  periods,  for  example,  from  about 
0lU5  to  0210  UT  and  from  about  1000  to  1630  UT  in  Figure  1.  The  termin¬ 
ation  of  the  radiation  often  appears  to  result  from  a  decrease  in  the 
bandwidth,  with  the  upper  and  lower  cutoffs  merging  as  the  radiation  dis¬ 
appears. 

If  the  radio  emission  is  generated  by  Doppler-shifted  cyclotron 
radiation,  as  in  the  models  of  Melrose  [1976!  and  Wu  and  Lee  [ 19791 »  a 
simple  explanation  can  be  advanced  for  the  upper  and  lower  cutoff  fre¬ 
quencies.  Before  discussing  the  origin  of  these  cutoffs,  it  is  useful  to 
first  review  the  essential  features  of  the  cyclotron  resonance  mechanism. 
From  very  general  considerations  the  growth  rate  for  the  right-hand 
extraordinary  mode  can  be  written 
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where  and  Nq  are  trie  hot  ana  cola  electron  densities,  to  and  w  are 
the  electron  plasma  frequency  and  gy rofrequencies ,  and  vj  and  Vj_  are  the 
parallel  and  perpendicular  velocities  relative  to  the  magnetic  field. 

The  integral  must  ue  evaluated  over  a  contour  C  in  velocity  space  defined 
by  the  cyclotron  resonance  condition. 


where  n  is  the  index  of  refraction,  6  is  the  wave  normal  angle  with 
respect  to  the  magnetic  field  and  u>  is  tne  wave  frequency.  As  pointed 
out  by  Wu  and  Lee  l 1979 1 ,  it  is  important  to  keep  the  relativistic  term 
in  Equation  2  even  for  velocities  much  less  than  the  speed  of  light. 

Two  free-energy  sources  can  be  identified  in  Equation  1,  one  corre¬ 
sponding  to  regions  of  positive  3F/3vp  and  the  other  corresponding  to 
regions  of  positive  3  K/  3  vjj  •  The  condition  3  F/  3  Vj_  >  0  is  characteristic  of 
a  ioss-cone  distribution  function,  and  3  F/ 3  v n >  0  is  characteristic  of  a 
beam- type  distribution.  A  loss-cone  distrioution  was  the  basis  of  Wu  and 
Lee's  i lyT 9 1  theory  of  auroral  kilometrie  radiation  ana  a  beam-type  dis¬ 
tribution  was  the  basis  of  Melrose's  liy'fb)  theory.  These  two  types  of 
distribution  functions  are  illustrated  in  Figure  2.  In  both  cases  tne 
presence  of  a  field-aligned  electrostatic  potential  gradient  plays  a 
crucial  r ol<--  in  establishing  the  essential  features  of  the  distribution 
function.  As  discussed  by  Chiu  and  Schulz  [1978],  in  the  presence  of  a 
parallel  electric  field  the  velocity  distribution  is  divided  into  five 
regions  by  a  loss-cone  boundary,  shown  by  the  solid  hyperbola-shaped  lines 
in  Figure  2,  and  an  electrostatic  acceleration  boundary,  shown  by  the 
dashed  elliptical  curve.  The  hyperbola-shaped  loss  cone  arises  because  of 


the  acceleration  caused  by  the  parallel  electric  field  between  the  space¬ 
craft  and  the  atmosphere.  The  elliptical  electrostatic  acceleration 
boundary  represents  the  boundary  for  magnetospheric  electrons  accelerated 
by  parallel  electric  fields  above  the  spacecraft.  The  shaded  region 
bounded  by  the  hyperbola  and  the  ellipse  consists  of  particles  trapped  in 
the  potential  well  formed  by  the  mirror  force  below  the  satellite  and  the 
electrostatic  potential  barrier  above  the  satellite.  For  the  loss-cone 
instability  it  is  sufficient  for  the  moment  to  point  out  that  a  large 
3F/3vp  is  expected  along  the  loss-cone  boundary.  In  the  case  of  the  beam- 
type  instability  a  large  3 F / 3 v u  arises  because  of  the  acceleration  of  low 
energy  magnetospheric  electrons  by  parallel  electric  fields  above  the 
spacecraft. 

Inspection  of  distribution  functions  measured  at  various  points  along 
the  auroral  field  lines  shows  evidence  of  both  the  beam  and  loss-cone 
free-energy  sources.  Figure  3,  for  example,  shows  a  typical  electron 
velocity  distribution  function  from  S3-3  [Mizera  and  Fennell,  1977]  which 
has  regions  of  large  3F/3vn  and  3F/3vp.  The  region  of  large  3F/3v||  occurs 
for  downward  moving  electrons  Just  inside  the  elliptical  electron  acceler¬ 
ation  boundary  (dashed  line)  and  inside  the  loss  cone.  The  distribution 
function  in  this  region  is  somewhat  different  than  the  bimaxwellian  used 
by  Melrose  for  the  electrostatically  accelerated  electron  beam,  but  the 
basic  free-energy  source  is  the  same.  The  region  of  large  3F/3vn  occurs 
at  a  velocity  of  approximately  vj  =  /  2e/b  s/m' ,  where  4S  is  the  electro¬ 
static  potential  at  the  spacecraft  measured  with  respect  to  the  potential 
above  the  acceleration  region.  The  large  parallel  velocity  gradient  in 
this  region  arises  because  of  the  markedly  decreased  particle  population 
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inside  the  electrostatic  acceleration  and  loss-cone  boundaries.  This 
region  of  velocity  space  is  only  accessible  to  electrons  from  the  iono¬ 
sphere,  which  have  much  lower  energies  and  therefore  much  lower  phase- 
space  densities  than  the  magnetospheric  electrons.  As  can  be  seen  in 
Figure  3,  the  region  of  large  3F/3vj_  occurs  for  upward  moving  electrons 
just  inside  the  loss-cone  boundary.  The  perpendicular  velocity  gradi¬ 
ents  are  particularly  large  near  and  inside  the  elliptical  electron 
acceleration  boundary,  at  parallel  velocities  V||  <  /  2e$  s/m* .  Examin¬ 
ation  of  other  distribution  functions  shows  that  this  relationship 
appears  to  be  a  consistent  feature  of  the  S3-3  observations.  The  V||  <■ 
/2^rs  Jm  constraint  on  the  region  of  large  3F/3vq  apparently  arises 
because  of  the  increased  electron  intensities  in  the  loss  cone  at  higher 
energies,  probably  due  to  backscattered  beam  electrons,  which  tend  to 
have  velocities  greater  than  / 2e$s/m' .  At  the  present  time  it  is  not 
known  which  of  the  two  free-energy  sources  is  most  important.  A  deci¬ 
sion  on  this  question  awaits  a  detailed  evaluation  of  the  growbh  rate 
from  Equation  1  using  the  measured  distribution  function.  Simple  con¬ 
siderations,  however,  suggest  that  the  loss-cone  source  is  probably  more 
important  because  the  factors  v^/c  and  CosG  strongly  reduce  the  contri¬ 
bution  due  to  3F/3v||  in  Equation  1. 

Having  discussed  the  essential  features  of  the  electron  distribu¬ 
tion  function  required  for  instability  we  can  now  discuss  the  physical 
reasons  for  upper  and  lower  cutoffs  in  the  radio  emission  spectrum.  In 
order  to  have  an  instability  it  is  essential  that  the  cyclotron  reson¬ 


ance  velocity  given  by  Equation  2  be  at  sufficiently  low  velocities  to 
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be  in  the  region  of  the  distribution  function  where  the  velocity  gradients 
are  large.  At  low  velocities,  v  <<  c,  and  for  an  index  of  refraction  n  - 
1,  Wu  and  Lee  [ 19T9 ]  have  shown  that  the  cyclotron  resonance  condition  is 
a  circle  in  the  V||  ,  vj_  plane  with  a  center  at  vn  =  cCos6  and  a  radius  of 
vR  =  c[Cos0  -  2Aw/wg)]1/2,  where  Aw  =  w-wg.  A  representative  resonance 
circle  is  shown  in  Figure  2.  Instability  occurs  when  w  and  G  are  such 
that  the  resonance  circle  passes  through  a  sufficiently  large  region  of 
positive  3F/8v||  or  3F/3vj_,  thereby  giving  a  positive  value  for  the  inte¬ 
gral  in  Equation  1.  It  is  evident  that  a  low  resonance  velocity  implies 
propagation  angles  close  to  tt/2  and  wave  frequencies  close  to  the  elec¬ 
tron  gyrof requency.  From  Equation  2  one  can  show  that  the  minimum  reson¬ 
ance  velocity  occurs  at  V||pes(Min)  =  c(Aw/w).  The  minimum  resonance  velo¬ 
city  is  therefore  a  very  sensitive  function  of  the  wave  frequency, 
decreasing  rapidly  as  w  approaches  wg.  The  wave  frequency,  however,  is 
constrained  to  be  above  the  extraordinary  mode  cutoff  which  occurs  at 


The  minimum  resonance  energy  is  then  approximately 


WHRES(Min)  =  Imc2 


Aw 


R=0 


(4) 


where  AwR_0  =  wR=0  -  wg  =  -wg/2  +  Awg/2)2  +  wp2'. 

The  variation  of  the  minimum  resonance  energy  with  radial  distance 
along  an  auroral  field  line  is  illustrated  in  Figure  4.  As  can  be  seen 
from  Equation  3  the  extraordinary  mode  cutoff,  wr=q,  tends  to  be  at  either 


oip  or  u)g,  which  ever  is  larger.  When  Wp  >  Ug  the  resonant  energies  are 
very  large  because  a  very  large  velocity  is  required  to  Doppler  shift  the 
emission  frequency  above  the  extraordinary  mode  cutoff.  As  Up  decreases 
below  oig,  the  minimum  resonance  energy  decreases  rapidly,  to  approximately 
W||  Res  (Min )  ~  (l/2)mc2  (cjp/wg)4  when  ujp  <<  ug.  As  shown  in  Figure  h  the  con¬ 
dition  dip  <  oig  tends  to  occur  at  intermediate  altitudes  between  the  iono¬ 
sphere  and  the  outer  magnetosphere.  If  the  resonance  energy  drops  below 
the  characteristic  energy  e^g  of  the  free-energy  source  in  the  electron 
distribution  (where  the  velocity  gradients  are  large)  then  instability  can 
occur.  For  a  given  electrostatic  potential  distribution  along  the  field 
line  this  condition  implies  upper  and  lower  boundaries  for  the  unstable 
region,  for  example,  A-A'  and  B-B^  in  Figure  4.  These  boundaries  in  turn 
imply  upper  and  lower  cutoffs  for  the  radio  emission  spectrum.  In  prac¬ 
tice,  of  course,  these  boundaries  are  not  as  precisely  defined  as  indicated 
in  Figure  U  because  the  growth  rate  given  by  Equation  1  involves  an  inte¬ 
gral  over  the  resonance  velocity.  This  model  does,  however,  provide  us 
with  a  qualitative  guide  to  analyze  how  the  bandwidth  of  the  emitted  radi¬ 
ation  depends  on  the  parameters  of  the  system. 

We  now  return  to  the  interpretation  of  the  observed  cutoffs  of  the 
auroral  kilometric  radiation.  It  is  evident  from  Figure  that  two  basic 
factors  control  the  bandwidth  of  the  emitted  radiation:  (l)  the  plasma 
density  distribution  via  the  dependence  of  W||pes(Min)  on  w p/wg,  and  (2) 
the  detailed  profile  of  the  electrostatic  potential  distribution  4>s.  Two 
extreme  situations  can  be  identified.  If  the  main  electron  acceleration 
occurs  at  very  high  altitudes,  as  in  case  B-B'  ,  then  the  bandwidth  is 
controlled  almost  completely  by  the  points  where  the  plasma  frequency 
profile  crosses  the  electron  gyrof requency  profile.  In  this  case,  the 


cutoff  frequencies  are  almost  completely  determined  by  the  plasma  density 
distribution.  On  the  other  hand  if  the  main  electron  acceleration  occurs 
at  lower  altitudes,  as  in  case  k-A' ,  then  the  lower  cutoff  frequency  is 
controlled  by  the  electron  gyrof requency  at  the  upper  boundary  of  the 
acceleration  region.  The  upper  cutoff  is  determined  by  a  combination  of 
factors.  If  the  potential  intersects  on  the  steep  part  of  the  W||ges 
curve  then  the  upper  cutoff  is  essentially  determined  by  the  plasma  den¬ 
sity  profile  at  the  point  where  o)p  =  wg.  If  the  potential  is  lower  then 
the  upper  cutoff  frequency  depends  on  the  details  of  the  potential  pro¬ 
file,  tending  to  decrease  as  $g  decreases. 

Not  much  is  known  about  the  variability  of  the  plasma  density  pro¬ 
file  at  high  altitudes.  The  main  information  comes  from  the  S3-3  satel¬ 
lite  at  about  2.0  Rg  and  various  eccentric  orbiting  spacecraft  at  some¬ 
what  greater  radial  distances.  Present  indications  are  that  u>p  <<  wg  at 
2.0  Rg  (R.  Torbert,  personal  communication),  more  or  less  in  agreement 
with  the  model  in  Figure  4.  At  larger  radial  distances,  >  5  Re*  the 
plasma  density  is  quite  low,  Ne  <  10  cm-3 ,  which  implies  that  the  plasma 
frequency  profile  probably  stays  below  the  electron  gyrof requency  out  to 
about  5  Rg ,  as  in  Figure  4.  Since  the  low  frequency  cutoff  of  the  kilo- 
metric  radiation  usually  corresponds  to  radial  distances  less  than  5  Re* 
it  is  tempting  to  assume  that  variations  in  the  low  frequency  cutoff  are 
caused  by  variations  in  the  upper  altitude  limit  of  the  acceleration 
region,  as  in  case  A-A'  of  Figure  4.  Unfortunately,  so  little  is  known 
about  the  density  distribution  along  the  auroral  field  lines  one  cannot 
be  certain  whether  this  assumption  is  valid  or  not.  If  it  is,  then  the 
low  frequency  cutoff  gives  a  direct  indication  of  the  upper  altitude 
limit  of  the  acceleration  region. 
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A  probability  distribution  for  the  upper  and  lower  cutoff  frequen¬ 
cies  obtained  from  a  portion  of  the  event  in  Figure  1  is  shown  in  Figure 
5.  A  radial  distance  scale  derived  from  the  electron  gyrofrequency, 
which  is  always  very  close  to  the  emission  frequency,  is  shown  at  the 
bottom  of  this  illustration.  If  the  plasma  frequency  remains  suffi¬ 
ciently  low  at  high  altitudes,  as  in  case  A-A' ,  these  low  frequency  cut¬ 
offs  imply  an  upper  limit  to  the  electron  acceleration  region  of  about 
2.0  to  3*5  Re*  Alternatively,  these  low  frequency  cutoffs  may  be  repre¬ 
sentative  of  the  radial  distance  at  which  the  plasma  frequency  crosses 
the  gyrofrequency,  as  in  case  B-B' .  The  interpretation  of  the  high  fre¬ 
quency  cutoff  is  more  certain.  At  frequencies  above  300  kHz  the  high 
frequency  cutoff  is  almost  certainly  controlled  by  the  crossover  between 
the  plasma  frequency  and  gyrofrequency  profiles  which  occurs  at  about 
1.5  Re*  At  lower  frequencies  this  cutoff  becomes  more  dependent  on  the 
electrostatic  potential  near  the  crossover  point,  with  a  decreasing 
upper  cutoff  frequency  indicating  a  smaller  accelerating  potential. 

As  discussed  earlier,  the  radiation  sometimes  disappears  completely 
during  certain  periods.  At  the  beginning  and  end  of  these  periods,  the 
upper  and  lower  cutoffs  appear  to  merge  as  though  the  bandwidth  of  the 
emission  decreased  to  zero  at  these  points.  Again  two  extreme  interpre¬ 
tations  are  possible.  The  termination  and  onset  of  the  emission  could 
be  due  to  changes  in  the  plasma  density  profile  or  to  changes  in  the 
electrostatic  potential  distribution.  If  changes  in  the  plasma  density 
are  the  controlling  effect  then  the  emissions  would  terminate  completely 
as  soon  as  the  plasma  frequency  exceeds  the  gyrofrequency  at  all  points 
along  the  magnetic  field  line.  If  such  a  density  control  does  occur  it 


seems  almost  certain  that  it  would  be  imposed  by  density  variations  at 
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high  altitudes  rather  than  by  variations  in  the  ionosphere.  Since 
rather  high  densities,  >  102cm-3  ,  would  be  required  in  the  outer  magnet¬ 
osphere  to  completely  quench  the  instability,  it  seems  unlikely  that 
changes  in  the  density  profile  could  play  the  dominant  role  in  switching 
the  radiation  on  and  off.  Temporal  variations  in  the  accelerating 
potential  appear  to  be  a  more  likely  alternative.  This  viewpoint  is 
supported  by  the  fact  that  the  kilometric  radiation  tends  to  disappear 
during  periods  of  low  magnetic  activity  [Voots  et  al. ,  19771  when  field- 
aligned  currents  and  the  accelerating  potential  are  expected  to  be  the 
smallest.  As  can  be  seen  from  Figure  H,  the  basic  model  predicts  that 
the  upper  and  lower  cutoffs,  A  and  A' ,  approach  each  other  as  the 
accelerating  potential  decreases,  thereby  shrinking  the  bandwidth  to 
zero  at  the  termination  and  onset  of  the  emission,  as  is  observed. 


III.  FINE  STRUCTURE 


Recently  wideband  measurements  from  the  ISEE  1  and  2  spacecraft  have 
provided  very  high  resolution  spectrum  measurements  of  auroral  kiloraetric 
radiation  [Gurnett  et  al. ,  19791*  The  spectrums  obtained  are  remarkable 
because  they  show  that  the  radiation  consists  of  many  discrete  narrowband 
emissions,  rather  than  a  continuous  broadband  emission  as  suggested  by 
the  lower  resolution  measurements  in  Figures  1.  Two  representative  high- 
resolution  frequency-time  spectrograms  of  auroral  kilometric  radiation 
are  shown  in  Figure  6.  These  spectrograms  cover  a  frequency  range  of  UO 
kHz,  from  125  to  165  kHz,  and  a  time  interval  of  about  20  minutes.  As 
can  be  seen,  the  spectrum  is  extremely  complex.  The  most  striking  char¬ 
acteristic  is  the  occurrence  of  many  narrowband  emissions  with  a  band¬ 
width  of  1  kHz  or  less.  The  entire  auroral  kilometric  radiation  spectrum 
appears  to  be  made  up  of  many  discrete  narrowband  tones.  Often  the  center 
frequency  of  the  individual  emissions  varies  in  a  systematic  manner, 
sweeping  either  upward  or  downward  across  the  spectrum.  The  spectrum  in 
Figure  6,  for  example,  is  dominated  by  emissions  drifting  upward  in  fre¬ 
quency  with  increasing  time.  Comparable  spectrums  can  also  be  found  in 
which  the  emissions  are  sweeping  downward  in  frequency.  Careful  examina¬ 
tion  shows  that  the  drifting  features  are  not  always  continuous  tones,  but 
sometimes  consist  of  short  bursts  occurring  on  time  scales  of  only  a  few 
seconds . 

A  contrasting  example  in  which  the  drifting  features  are  nearly 
absent  is  illustrated  in  Figure  7,  which  shows  a  UO-kHz  frequency  range, 
from  62.5  to  102.5  kHz,  for  two  successive  1-minute  intervals.  The 


spectrum  in  this  case  consists  of  a  very  large  number  of  brief  narrowband 
bursts.  The  duration  and  frequency  range  of  the  individual  bursts  is 
very  small,  typically  1  sec  duration  and  3  kHz  frequency  range,  or  less. 
The  spectral  shapes  of  the  individual  bursts  vary  over  a  wide  range, 
sometimes  sweeping  upwards  or  downwards  in  frequency  and  forming  various 
curved  features  on  the  frequency-time  diagram. 

The  extremely  complex  spectral  structure  of  the  auroral  kilometric 
radiation  presents  some  difficult  interpretational  questions.  The  occur¬ 
rence  of  many  discrete  narrowband  emissions,  as  in  Figure  7,  is  strongly 
reminiscent  of  chorus  and  other  types  of  discrete  whistler-mode  emissions 
observed  in  the  earth's  magnetosphere  [Helliwell,  1965]  .  The  close 
similarity  of  the  dynamic  spectral  characteristics  may  be  more  than 
purely  coincidental,  since  both  the  auroral  kilometric  and  whistler-mode 
instability  mechanisms  may  be  of  essentially  the  same  type.  If  the  kilo¬ 
metric  radiation  is  produced  by  a  loss  cone  distribution  then  both  insta¬ 
bilities  involve  a  cyclotron  resonance  interaction  with  an  anisotropic 
electron  distribution.  The  detailed  resonance  condition  and  index  of 
refraction  are,  of  course,  quite  different,  but  otherwise  the  physical 
process  is  quite  similar.  For  whistler-mode  emissions  it  is  generally 
agreed  that  the  dynamic  evolution  of  the  spectrum  is  a  highly  nonlinear 
process  involving  electrons  trapped  by  the  rotating  wave  field  llnan  et 
al. ,  19781-  It  is  possible  that  similar  processes  may  be  able  to  account 
for  the  discrete  structure  of  the  auroral  kilometric  radiation.  It  is 
interesting  to  note  that  if  the  generation  of  auroral  kilometric  radi¬ 
ation  is  essentially  similar  to  the  whistler-mode  instability,  then  the 
kilometric  radiation  may  play  a  significant  role  in  the  pitch-angle 
scattering  of  auroral  electrons,  similar  to  the  pitch-angle  scattering  of 
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radiation  belt  electrons  by  whistler-mode  emissions  in  the  inner  regions 
of  the  magnetosphere.  Because  the  auroral  kilometric  radiation  is 
generated  in  the  right-hand  polarized  extraordinary  mode  the  radiation 
carries  away  right-hand  angular  momentum,  thereby  reducing  the  pitch 
angle  of  the  interacting  electrons,  very  similar  to  the  situation  with 
whistler-mode  interactions.  Because  the  radiation  intensities  in  the 
source  are  undoubtedly  very  large  { 10  to  100  mVm-1  ),  the  pitch-angle 
diffusion  rates  could  be  a  very  significant  factor  in  controlling  the 
equilibrium  electron  velocity  distribution  in  the  auroral  acceleration 
region.  If  such  scattering  occurs  it  should  be  characterized  by  time 
scales  comparable  to  the  time  scale  of  the  individual  bursts,  which  is 
only  a  few  seconds.  Rapid  fluctuations  in  inverted-V  electron  precipi¬ 
tation  fluxes  have  been  reported  by  various  investigators,  including  for 
example  Lin  and  Hoffman  [1979]  ,  which  could  possibly  be  related  to  the 
fine  structure  of  the  kilometric  radiation. 

Probably  the  most  important  single  feature  which  must  be  explained 
in  the  fine  structure  of  the  auroral  kilometric  radiation  is  the  organi¬ 
zation  of  the  discrete  bursts  into  upward  and  downward  drifting  bands  of 
the  type  illustrated  in  Figure  6.  The  occurrence  of  quasi-monochromatic 
emissions  sweeping  through  a  large  range  of  frequencies  is  a  common  fea¬ 
ture  of  solar  radio  bursts  [Kundu,  1965]  and  Jovian  decametric  radio 
emissions  [Warwick,  I96T ] -  These  drifting  features  are  usually  inter¬ 
preted  as  being  due  to  the  motion  of  the  emitting  particles,  or  a  propa¬ 
gating  disturbance,  through  a  plasma  which  has  a  spatial  gradient  in 
the  characteristic  emission  frequency.  This  model  is  used,  for  example, 
to  explain  the  characteristic  frequency  drifts  of  type  III  and  type  IV 
solar  radio  bursts.  In  the  case  of  the  auroral  kilometric  radiation  the 
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characteristic  emission  frequency  is  the  electron  gyrofrequency.  Since  the 
resonance  energy  is  very  sensitive  to  the  frequency  difference.  Aw,  between 
the  emission  frequency  and  the  electron  gyrofrequency,  it  is  difficult  to 
see  how  one  can  account  for  the  large  frequency  range  of  the  drifting  emis¬ 
sions  without  involving  a  motion  of  the  source. 

Given  that  the  emission  takes  place  near  the  gyrofrequency,  it  is 
relatively  easy  to  estimate  the  radial  component  of  the  source  velocity, 
dR/dt,  from  the  observed  frequency  drift  rate,  df/dt.  For  motion  along  a 
dipole  field  line  at  high  latitudes,  the  velocity  of  the  source  is  approx¬ 
imately 
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where  fgo  =  1.7  MHz  is  the  electron  gyrofrequency  at  the  surface  of  the 
earth  in  the  auroral  region.  Using  this  equation  we  have  estimated  the 
source  velocity  using  a  random  sample  of  the  drifting  emissions  observed  in 
the  frequency  range  from  125  to  500  kHz.  The  distribution  of  source  velo¬ 
cities  is  illustrated  in  Figure  8.  As  can  be  seen,  the  source  velocities 
vary  over  a  large  range,  from  about  3  to  300  km/sec,  with  a  median  value  of 
about  30  km/sec.  Downward  source  motions  are  much  more  common  than  upward 
source  motions. 

To  aid  in  the  interpretation  of  these  apparent  source  velocities  cor¬ 
responding  electron  and  proton  energy  scales  are  shown  at  the  top  of 
Figure  8.  As  can  be  seen,  the  observed  source  velocities  are  much  smaller 
than  the  typical  electron  and  proton  energies,  100  eV  to  10  keV,  observed 
along  the  auroral  field  lines.  The  observed  drift  rates  cannot,  therefore, 
be  identified  with  the  actual  velocities  of  the  auroral  particles.  In 
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considering  other  characteristic  velocities  which  could  be  involved,  two 
possibilities  immediately  come  to  mind:  the  Alfven  speed,  =  B  //  y  Qp m',  and 
the  ion-acoustic  speed,  Vq  =  / kTe/mj’.  Because  the  plasma  density  is  quite 
low  in  the  resonance  interaction  region  (wp  <  Log)  the  Alfven  speed  is  too 
large,  >  2  x  103  km/sec,  to  be  considered.  The  ion-acoustic  speed  is 
mainly  determined  by  the  electron  temperature,  Te,  and  the  ion  mass,  m^. 
Using  a  temperature  of  104  °K  and  assuming  that  0+  is  the  dominant  ion 
species,  the  ion-acoustic  speed  in  the  ionosphere  is  estimated  to  be  about 
Vc  *  2.3  km/sec.  For  H+  the  ion-acoustic  speed  is  about  Vc  =  9*2  km/sec. 

At  higher  altitudes  the  ion-acoustic  speed  tends  to  increase.  An  upper 
limit  is  given  by  the  ion-acoustic  speed  in  the  plasma  sheet,  which  is 
approximately  Vq  =  200  km/sec.  This  range  of  ion-acoustic  speeds  is  indi¬ 
cated  in  Figure  8.  As  can  be  seen,  the  observed  source  velocities  are  in 
approximately  the  same  range  as  the  estimated  ion-acoustic  speeds.  This 
agreement  suggests  that  the  drifting  tones  observed  in  the  auroral  kilo- 
metric  radiation  may  be  caused  by  shock-like  disturbances  propagating  along 
the  auroral  field  lines  at  the  ion-acoustic  speed.  The  shock-like  char¬ 
acter  of  the  disturbance  is  indicated  by  the  narrow  bandwidth  of  the  drift¬ 
ing  emissions,  which  suggests  a  spatial  scale  for  the  disturbance  of  about 
AR  =  50  km.  This  mechanism  would  then  be  somewhat  similar  to  the  genera¬ 
tion  of  type  IV  solar  radio  bursts  which  are  produced  by  shock  waves  propa¬ 
gating  through  the  solar  corona  [Kundu,  I965I •  It  has  long  been  suspected 
from  computer  simulations  [Hubbard  and  Joyce,  1979!  and  laboratory  experi¬ 
ments  [Carlqvist  and  Bostrom,  19T0]  that  electrostatic  structures  in  the 
auroral  acceleration  region  may  be  highly  turbulent,  possibly  consisting  of 
many  short  duration  transient  disturbances  propagating  at  speeds  near  the 


ion-acoustic  speed.  The  exact  mechanism  by  which  these  disturbances  could 


control  the  spectrum  of  the  auroral  kilometric  radiation  is  not  known, 
although  numerous  possibilities  exist.  If  the  potential  perturbations 
are  sufficiently  large,  then  the  disturbances  may  significantly  modify 


the  local  electron  velocity  distribution,  either  by  changing  the  loss- 
cone  boundary,  via  the  dependence  on  $s,  or  by  changing  the  reflection 
condition  for  the  electrostatically  trapped  electron  distribution.  The 
changes  could  in  turn  affect  the  growth  rate  and  spectrum  of  the  radio 
emission,  thereby  producing  drifting  spectral  features  which  appear  to 
move  with  the  disturbance. 


IV.  CONCLUSION 


In  this  paper  we  have  discussed  the  principal  features  of  the  spec¬ 
trum  of  auroral  kilometric  radiation.  Using  current  theories  we  have 
tried  to  relate  these  features  to  processes  occurring  in  the  auroral 
acceleration  region.  At  the  present  time,  it  appears  that  the  upper  and 
lower  cutoff  frequencies  of  the  auroral  kilometric  radiation  can  be 
readily  understood  if  the  radiation  is  produced  by  a  Doppler-shifted 
cyclotron  emission  driven  by  either  a  loss-cone  or  beam  distribution. 

The  interpretation  of  the  detailed  variations  of  the  upper  and  lower 
frequency  cutoffs  is  complicated  by  uncertainties  about  the  relative 
importance  of  variations  in  the  plasma  density  profile  and  variations  in 
the  electrostatic  potential  distribution  along  the  auroral  field  lines. 
If  the  electron  plasma  frequency  is  sufficiently  low,  £  20  kHz,  at  high 
altitudes  in  the  magnetosphere,  then  the  lower  frequency  cutoff  of  the 
radio  emission  spectrum  gives  the  upper  limit  to  the  auroral  acceler¬ 
ation  region.  Typically,  the  upper  limit  of  the  acceleration  region  is 
located  at  radial  distances  from  about  2.0  to  3-5  Rg*  During  some 
intervals  the  kilometric  radiation  completely  disappears.  These  inter¬ 
vals  are  believed  to  correspond  to  either  periods  when  the  plasma  den¬ 
sity  is  unusually  high,  at  all  points  along  the  auroral  field 

lines,  or  (more  likely)  to  periods  of  very  small  accelerating  poten- 


21 


High  resolution  frequency-time  spectrums  of  auroral  kilometric 
radiation  show  that  the  radiation  consists  of  many  discrete  narrowband 
emissions,  some  of  which  drift  upwards  or  downwards  over  a  large  fre¬ 
quency  range.  The  existence  of  this  very  complex  fine  structure  indi¬ 
cates  the  presence  of  corresponding  complex  spatial  and  temporal  struc¬ 
ture  in  the  auroral  acceleration  region.  The  drifting  narrowband  emis¬ 
sions  indicate  the  presence  of  very  compact  source  regions,  with  radial 
scale  sizes  of  only  50  km,  which  move  upward  or  downward  along  the  mag¬ 
netic  field  line  with  velocities  ranging  from  3  to  300  km/sec.  These 
source  velocities  are  comparable  to  the  ion-acoustic  speed  in  the 
auroral  ionosphere,  which  suggests  the  existence  of  shock-like  distur¬ 
bances  propagating  along  the  auroral  field  lines.  These  disturbances 
may  possibly  originate  from  unstable  fluctuations  in  the  electrostatic 
potential  distribution,  possibly  forming  transient  double-layers  or 
shock-like  structures  propagating  along  the  magnetic  field  line.  The 
fine  structure  of  the  auroral  kilometric  radiation  also  has  features 
which  are  remarkably  similar  to  discrete  whistler-mode  emissions,  sug¬ 
gesting  that  somewhat  similar  nonlinear  processes  may  be  involved  in 
both  types  of  radio  emissions. 
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Figure  1 


Figure  2 


Figure  3 


FIGURE  CAPTIONS 


A  representative  spectrum  of  auroral  kilometric 
radiation.  Note  the  upper  and  lower  cutoff 
frequencies  of  the  spectrum  and  the  occurrence  of 
intervals  when  the  radiation  disappears  completely. 

Idealized  velocity  distributions  illustrating  the 
loss-cone  and  beam-type  free-energy  sources. 

Cyclotron  resonance  interactions  occur  along  a  circle 
in  velocity  space  with  a  center  at  vj  =  cCos0  and  a 
radius  v^  *  clCosB  -  2Au>/wgi  l'2. 

A  velocity  distribution  function  obtained  by  S3-3  in 
the  auroral  acceleration  region  [Mizera  and  Fennel, 
1977]*  Note  the  regions  of  large  positive  3F/3vq  near 
the  loss  cone  boundary  and  large  positive  3F/3v|  just 
inside  the  elliptical  electrostatic  acceleration 
boundary . 
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An  ionospheric  model  illustrating  the  control  of  the 
upper  and  lower  cutoffs  of  the  kilometric  radiation  by 
the  electron  density  profile  and  the  electrostatic 
potential  distribution.  Instability  only  occurs  if 
the  minimum  resonance  energy,  W||pgs(Min)>  is  below  the 
electrostatic  potential  energy  e4>s. 

The  distribution  of  upper  and  lower  cutoff  frequencies 
for  Figure  1,  and  the  corresponding  radial  distance, 
assuming  that  the  emission  occurs  very  close  to  the 
electron  gyrofrequency. 

High  resolution  spectrograms  showing  the  occurrence  of 
many  narrowband  emissions  drifting  upward  in 
frequency. 

High  resolution  spectrograms  showing  the  occurrence  of 
many  discrete  emissions  in  the  auroral  kilometric 
radiation  spectrum.  These  emissions  have  many 
features  similar  to  discrete  whistler-mode  emissions. 

The  distribution  of  drift  velocities  for  narrowband 
drifting  features  similar  to  those  illustrated  in 
Figure  7. 
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